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ABSTRACT 

For varying values of tan (3, we study the effect of CP violating phases from 
the soft supersymmetry breaking terms in string-inspired models on the 
relic abundance and detection rates of the lightest neutralino (LSP). We 
find that the phases have no significant effect on the LSP relic abundance 
but can have a substantial impact on the detection rates. 
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1. Introduction 

CP violation is an important test for physics beyond the standard model. In the 
standard model, only one CP violating phase exists in the Kobayashi-Maskawa matrix. 
However, in the supersymmetric standard model there are many complex parameters, in 
addition to Yukawa couplings, which lead to new sources of CP violation. These include 
the mass coefficient fi of the bilinear term involving the two Higgs doublets, the SU{3), 
SU{2) and U{1) gaugino masses M3, M2 and Mi, and the parameters Af and B which 
respectively are the coefficients of the supersymmetric breaking trilinear and bilinear 
couplings. ( the subscript / denotes the flavor index.) 

In minimal supersymmetric standard model (MSSM), only two of these phases are 
physical. Through appropriate field redefinitions we end up with the phase of fi {(p^) and 
the phase of A {(Pa) as the physical phases which cannot be rotated away. The phase 
of B is fixed by the condition that Bfi is real. It is known that, unless these phases are 
sufficiently small, their contributions to the neutron electric dipole moment (EDM) are 
larger than the experimental limit 1.1 x 10~^^ e.cm. Recently, the effect of these phases on 
the EDM of the neutron was examined in a model with dilaton-dominated supersymmetry 
(SUSY) breaking, taking into account the cancellation mechanism between the different 
contributions. It was shown that for a wide region of the parameter space the phase of n 
is constrainted to be of order 10~^, while the phase of A is strongly correlated with that 
of fi in order not to violate the bound on the neutron EDM. 

The effect of SUSY CP violating phases on the relic density of the LSP has been 
considered for the MSSM case in Ref. |]12|, and for the supersymmetric standard model 



coupled to N=l supergravity in Ref. |jT3[. It was shown in Ref. [|T2| that the upper 
bound on the LSP (from Qh"^ < 0.25) is relaxed from 250 GeV to 650 GeV. The effect 
of CP violation on the direct detection rates of the LSP in MSSM is also considered in 



Ref. [^, 0. We argue that such a large upper bound on the LSP is not possible in the 
model we consider here. We show that in case of the bino-like LSP the chance of the CP 
phases to have a significant effect is very small. The impact of the phases on the direct 
and indirect detection rates is an important issue and we present some details here. 

The paper is organized as follows. In section 2 we discuss the effect of the CP phases 
on the LSP mass and purity within the string inspired model considered in Ref. 0. In 
section 3 we compute the relic abundance of the LSP for low and intermediate values of 
tan f3. We find that the CP phases have almost no effect on the LSP relic density, so that 
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the upper bound on the LSP mass obtained in Ref. |[TT|, ^ remains unchanged. In section 
4 we discuss the large tan/5(~ rrit/mf,) case and again find that there is no significant 
effect of CP phases on the LSP rehc density. In section 5 we show that CP phases can 
have a substantial effect on the LSP detection rates. Our conclusions are given in section 
6 



2. String inspired model 

We will consider the string inspired model which has been recently studied in Ref. . In 
this model, the dilaton S and overall modulus field T both contribute to SUSY breaking. 
The soft scalar masses and the gaugino masses Ma are given as 

mf = 7713/2(1 + cos^ 6^), (1) 
Ma = v^m3/2sin^e-'°^ (2) 

where m3/2 is the gravitino mass, rii is the modular weight of the chiral multiplet, and 
sin 6' defines the ratio between the F-terms of S and T, (For example, the limit sin 6' — > 1 
corresponds to a dilaton-dominant SUSY breaking). The phase as originates from the 
F-term of S. 

The A-terms can be written as 

Ajk = -\/3m3/2 sin ^e"'°^ - 7773/2 cos 6^(3 + rij + r7j + r7fc)e"*°^, (3) 

where n,, rij and are the modular weights of the fields that are coupled by this A-term. 
One needs a correction term in eq (^ when the corresponding Yukawa coupling depends 
on moduli fields. However, the T-dependent Yukawa coupling includes a suppression 
factor 1^, and so we ignore it. Finally, the phase originates from the F-term of T. 

The magnitude of the soft SUSY breaking term BfiHiH2 depends on the way one 
generates a 'natural' /x-term. Here we take fi and B as free parameters and we will fix 
them by requiring successful electroweak (EW) symmetry breaking. 

As stated earlier, the gaugino masses as well as A-terms and 5-term are, in general, 
complex. We have the freedom to rotate Ma and Aij^ at the same time 0. Here we use 
the basis in which Ma is real. Similarly, we can rotate the phase of B so that Bfi itself 
is real. In other words, (pB = = arg(i?M*). In this basis, A-terms contain a single 
phase, {a A = ax — o^s)- 
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As shown in eqs.([I|-0), the values of the soft SUSY breaking parameters at string scale 
depend on the modular weights of the matter states. The modular weights of the matter 
fields Hi are normally negative integers. Following the approach of Ref. [0 the 'natural' 
values of modular weights for matter fields ( in case of orbifolds) are -1,-2,-3 and -4. It 
was shown in Ref. |^ that the following modular weights for quark and lepton superfields 
is favorable for EW breaking 

nQ = nu = uh^ = -1, 

and 

ud = ul = He = = -2. 
Under this assumption we have 

At = Ab = -V3m3/2 sme + mg/a cos^e"*"^, (4) 

and 

Ar = - ^3^3/2 sin 6^ 2 mg/a cos 06^'°''' , (5) 



Given the boundary conditions in eqs. (0-0) at the compactification scale, we determine 
the evolution of the couplings and the mass parameters according to their one loop renor- 
malization group equation in order to estimate the mass spectrum of the SUSY particles 
at the weak scale. The radiative EW symmetry breaking imposes the following conditions 
on the renormalized quantities: 

rriH, + + > 2S/i, (6) 
iml^+fi'){ml^+fi')<{BfiY, (7) 

^ - tan2/3-l 2"' 

and 

-2Ba 

sin 2/? = , ^ , , (9 

m2^^ +m2^^ +2/i2' ^ ^ 

where tan /? = {H2) / {H^ is the ratio of the two Higgs VEVs that gives masses to the up 
and down type quarks, and m'jj_^, mfj^ are the two soft Higgs square masses at the EW 
scale. Using the above equations we can determine and B in terms of ms/2, d and ua- 
The phase of 0^ remains undetermined. 
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Since we are interested in investigating the effect of the supersymmetric phases and 
(f)^ on the rehc density of the LSP and its direct and indirect detection rates, we first 
study the allowed regions of these phases and later impose the constraints (derived in 
Ref 0) from the experimental bounds on the electric dipole moments. 

The neutralinos Xi, = 2, 3, 4) are the physical (mass) superpositions of the Higgsi- 
nos H^, H2 and the two neutral gaugino (bino) and (wino). The neutralino mass 
matrix is given by 

/ Ml -MzcoslSsinOw Mz sin ^ sin % \ 

M2 Mz cos 13 cos 9w —Mzsin(5cos9w 

viN = , , (10) 

—Mz cos f3 sinOw Mz cos P cos 9w ^e'''^' 

\ M^sin/3sin6'H/ - Mz sin (3 cos Ow /xe<^^ / 

where Mi and M2 now refer to 'low energy' quantities whose asymptotic values are given 
in equation (^. The lightest eigenstates Xi is a linear combination of the original fields: 

X? = NnB + Ni2W^ + + NuHl (11) 

where the unitary matrix Nij relates the Xi fields to the original ones. The entries of this 
matrix depend on 777.3/2, 6 and (p^. The dependence of the Xi (LSP) mass on (p^ is shown 
in figure 1 for 77^3/2 — lOOGeV, cos^ 9 ~ 1/2 and a a — 7r/2. 
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Figure 1. LSP mass as a function of the phase 0^ . 
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A useful parameter for describing the neutralino composition is the gaugino "purity" 
function 

/, = |iVup + |iVi2p (12) 

We plot this function versus 0^ in figure (2) which clearly shows that the LSP is essentially 
a pure bino. These two figures show that the neutralino mass and composition are only 
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Figure 2. LSP gaugino purity as a function of the phase (f)^ . 

slightly dependent on the supersymmetric phase (j)^. 

3. Relic Abundance Calculation for low and 

INTERMEDIATE VALUES OF tan/3 

In this section we compute the relic density of the LSP in the case of low tan (3 {i.e. 
tan/3 ~ 3), as well as for intermediate tan/5 values {i.e. tan/5 ~ 15). Using a standard 
method ||^ in which we expand the thermally averaged cross section (a^f ) as 

{aAv) = a + hv'^ + ... , (13) 

where v is the relative velocity, a is the s-wave contribution at zero relative velocity and 
h contains contributions from both the s and p waves, the relic abundance is given by . 

^x^' = = 2-82 X 10''Y^{mJGeV), (14) 

Pc/ iT' 
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where 

= 0.264 g'J' Mp (— + — ), (15) 

h is the well known Hubble parameter, 0.4 < h < 0.8, and pc ~ 2 x 10~^^/i^ is the critical 
density of the universe. The freeze-out temperature is given by 

Xf = In ^ , ' ^ '—^ (16) 

V9*xf 

119 



Here xp = m^/TF, Mp = 1.22 x 10 GeV is the Planck mass, and g^, (8 < < 10) is 



the effective number of relativistic degrees of freedom at Tp. Also c = 1/2 as explained 



in Ref flO 



Given that the LSP is bino-like, the annihilation is predominantly into leptons, with the 
other channels either closed or suppressed. It is worth noting that the squark exchanges 
are suppressed due to their large masses, as figure (3) shows. 
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Figure 3. The squarks (dashed-dotted) and sleptons (solid line) masses 
versus the LSP mass, for cos^6' ~ 1/2 and — 7r/2 . 



The annihilation process is dominated by the exchange of the right slepton. In fact, 
the masses of the right slepton are essentially independent of and 0^, unless there is 
a significant amount of slepton mixing. Here, the off-diagonal element of the matrices 
are much smaller than the diagonal elements M^^ and M^^. Furthermore, since the LSP 
is essentially a bino, it only slightly depends on the phase of fi as figure (1) confirms. 
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Therefore, we find that the constraint on the relic density: 0.1 < I^lsp ^ 0.9, with 
0.4 < h < 0.8 leads, as figure (4) shows, to the previously known upper bound on the 
LSP mass found in the case of vanishing SUSY phases [^, |T^, namely, < 250 GeV. 
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Figure 4. The relic abundance of the LSP versus its mass. 



This result is different from the one discussed in Ref. 12 where it was claimed that 



the CP violating phases have a significant effect in that the cosmological upper bound on 
the bino mass is increased from 250 GeV to 650 GeV. This enhancement can be traced 
to the assumptions proposed in that model, namely that all the scalar masses are equal 
and of order Mw at the weak scale. Also, the sfermion mxings were assumed to be large 
(with fi ~ TeV). It turns out that such assumptions may lead to unacceptable charge 
and color breaking, as explained in Ref. W3\. Moreover, they cannot be motivated from 
supergravity or superstring models. 

We have also considered the relic density for intermediate tan/3 {i.e. tan/? ~ 15). We 
find that there is no significant difference between this and the case of low tan (3. The 
upper bound on the LSP mass is still of order 250 GeV. 

4. SUSY CP PHASES WITH Large tan/? 

We now extend our study to the case when tan/5 is large (~ 50). In a large class of 
supersymmetric models with flavor U{1) symmetry, tan/5 ~ (^).e"', where e ^ 0.2 (~ 
Cabibbo angle) is a 'small' expansion parameter, and n = 0,1,2. (See for instance IIP 
and reference therein). 
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For tan/5 ~ it is known that the phenomenological aspects of these models are 
very different compared with the small tan jS case. In particular, radiative EW symmetry 
breaking is an important non-trivial issue. Non-universality, such as m|^^ > mf^^ ^^e 
Planck scale, is favored for a successful EW breaking with large tan f3. Further, non- 
universality of the squark and slepton masses can affect symmetry breaking as well as 
other phenomenological aspects. We have adopted this non-universality in our choice for 
the modular weights in section 2. 

In the large tan/3 case the Higgs potential has two characteristic features. It follows 
from the minimization conditions that 

ml ^ (17) 
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with 

~ + > 0. (19) 

Here, mf = m'jj_ + fi,i = 1,2 and m| = Bfi. A combination of eqs. ([TTI) and (0) gives the 
following constraint on the low energy parameters 

- > M|, (20) 

i.e — m'jj^ > M|. In order to have electroweak breaking in the large tan/3 case, the 
difference between the masses of the two Higgs fields should satisfy the above inequality. 

In our model we find that this inequality is indeed satisfied, and the EW symmetry is 
broken at the weak scale. Also, one of the stau leptons (f/j) has a 'small' mass of order 
O(IOO) GeV, and happens to be the lightest slepton. It therefore dominates the LSP 
annihilation process. This relaxes the upper bound on the LSP mass from 250 to 300 
GeV. Thus, even in the of large tan f3 case the effect of the supersymmetric phases are 
relatively small as figure (5) shows. This essentially follows because the diagonal elements 
of the stau mass matrices, in this model are larger than the off-diagonal ones, i.e., there 
is no large mixing, as well as from the fact that the LSP is bino like. 
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Figure 5. The relic abundance with non- vanishing phases (sohd hne) and 
vanishing phases (dashed hne) versus the LSP mass in case of tan/5 ~ 



5. CP PHASES AND DETECTION RATES OF THE LSP 

We have seen that the effects of CP violating phases on the neutralino relic density are 
very small. In this section we examine the effect of these phases on the event rates of relic 
neutralinos scattering off nuclei in terrestrial detectors. The direct detection experiments 
provide the most natural way of searching for the neutralino dark matter. Any large CP 
violating phases can affect the detection rate, as will see below. It is interesting to note 
that the measured event rate may shed light on the value of the super symmetric phases. 



The differential detection rate is given by 



dR 
dQ 



(21) 



where fi{v) is the distribution of speeds relative to the detector. The reduced mass is 



rrir 



where m^r is the mass of the nucleus, f„ 



S^)i/2, Q is the energj 



deposited in the detector, and is the density of the neutralino near the Earth, a is the 
elastic-scattering cross section of the LSP with a given nucleus. In general a has two con- 
tributions: a spin-dependent contribution arising from and q exchange diagrams, and 
a spin-independent (scalar) contribution due to the Higgs and squark exchange diagrams. 
For "^^Ge detector, where the total spin of ''^Ge is equal to zero, we have contributions 
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only from the scalar part. 



Ami 



a 



[ZU+{A-Z)Q\ 



(22) 



where Z is the nuclear charge, and A — Z is the number of neutrons. The expressions for 
/„ and /„, and their dependence on the SUSY phases can be found in Ref. O, llO]. The 



effect of the CP violating phases enter through the neutralino eigenvector components 
Nij, and also through the matrices that diagonalize the squark mass matrices. Finally, 
F[Q) in (|2l|) is the form factor. We use the standard parameterization |17 

3ji(gi?i) 



F{Q) 



qRi 



12 2 

-e ■5'?'^', 



(23) 



where the momentum transfer = 2mNQ, Ri = (R^ — 5s^)^/^ with R = 1.2fmA^^'^, and 
A is the mass number of "^^Ge. ji is the spherical Bessel function and s ~ Ifm. 

The ratio R of the event rate with non- vanishing CP violating phases to the event rate 
in the absence of these phases is presented in figure (6). The solid curve corresponds to 
the case = 0, while the dashed one corresponds to a a = 7r/2. 
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Figure 6. The ratio of the direct detection rates as function of 



From this figure, it is clear that, in the model we are considering, the CP violating 
phases can significantly affect the event rates of the direct detection of the LSP. The 
phase reduces the value of R, while the phase in the trilinear coupling increases 
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it. However, as explained in Ref. 0, 0^ is constrained from the electric dipole moment 
experimental limit to be < 10~^ 

For completeness, we also examine the effect of CP violating phases on the indirect 
detection rates of the LSP in the halo. The observation of energetic neutrinos from the 
annihilation of the LSP that accumulate in the sun or in the earth is a promising method 
for detecting them. The technique for detecting such energetic neutrinos is through obser- 
vation of upward going muons produced by charged current interactions of the neutrinos 
in the rock below the detector. The flux of such muons from neutralino annihilation in 
the sun is given by 

F = 2.9 X 10«m-V-Hanh^(t/r)pf /K)CK)(^)^(^^)^. (24) 
The neutralino-mass dependence of the capture rates is described by [§] 

/(mj = J2 /.0.g.(^x)i^.K) (^^;"'4p , (25) 

where the quantities (pi and fi describe the distribution of element i in the sun and they 
are listed in Ref. 0, the quantity Si{m^) = S{-^^) is the kinematics suppression factor 
for the capture of neutralino of mass from a nucleus of mass rriN^ 0, and Fi{m^) is 
the form factor suppression for the capture of a neutralino of mass by a nucleus i. 
Finally, the function C(^x) describes the energy spectrum from neutralino annihilation 
for a given mass. 



In Figure (7) we present the ratio of the muon fluxes resulting from captured neutralinos 
in the sun in the case of non- vanishing CP violating phases to that of vanishing 0^, for 

= 0.3GeV/cm^. We see that the predicted muon flux increases as the phase of A-term 
is increased. 

We can understand this important effect of CP violating phases on the detection rate 
as follows. The phases affect the neutralino eigenvector components Nij and the squark 
mass matrices. Consequently, they have a significant effect on the neutralino coupling 
to quarks. The spin independent contribution, as also shown in Ref. [|15| and |T^, is 
decreased by increasing the phase of fi, and goes in the other direction if the phase of 
A-term is increased. This leads to the same behaviour for the elastic scattering cross 
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Figure 7. The ratio of the indirect detection rates as function of (f)^ 

section, which translates this dependence on the phases of /j, and A to the detection rates 
as figures 6 and 7 confirm. 

6. Conclusions 

We have studied the impact of CP violating phases from soft SUSY breaking terms in 
string-inspired models on the LSP, its purity and its relic abundance density. For different 
values of tan/? (of order unity, intermediate and of order mt/m^), we found that these 
phases have no significant effect on the LSP relic density, so that the upper bound on the 
LSP mass is essentially unchanged. We have examined the effect of these phases on the 
direct and indirect detection rates. We found that increasing the value of the phase 0^ 
leads to a decrease the event rates, while the phase (pA of the trilinear coupling has the 
opposite effect. 
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